Luttinger liquid universality in the time evolution after an interaction quench 
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We provide strong evidence that the relaxation dynamics of one-dimensional, metallic Fermi 
systems resulting out of an abrupt amplitude change of the two-particle interaction has aspects which 
are universal in the Luttinger liquid sense: The leading long-time behavior of certain observables is 
described by universal functions of the equilibrium Luttinger liquid parameter and the renormahzed 
velocity. We analytically derive those functions for the Tomonaga-Luttinger model and verify our 
hypothesis of universality by considering spinless lattice fermions within the framework of the density 
matrix renormalization group. 
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The equilibrium low-energy physics of a large class of 
one-dimensional (Id), correlated, metallic Fermi systems 
is described by the Luttinger liquid (LL) phenomenology 
[H, 0] ■ The Tomonaga-Luttinger (TL) model is the effec- 
tive low-energy fixed point model of the LL universality 
class and thus plays the same role as the free Fermi gas in 
Fermi liquid theory. This universality relies on the renor- 
malization group (RG) irrelevance of contributions such 
as the momentum dependence of the two-particle interac- 
tion Q or the curvature of the single-particle dispersion 
[3| which are present in microscopic models but ignored 
in the TL model. For a model falling into the LL uni- 
versality class it is not necessary to explicitely compute 
thermodynamic observables and correlation functions if 
one is interested in the low-energy limit. One only needs 
to determine two numbers - the LL parameter K and 
the renormahzed velocity v of the excitations - which 
fully characterize the low-energy physics of a spinless LL 
(on which we focus). Those depend on the band struc- 
ture and filling as well as the amplitude and range of 
the two-particle interaction of the microscopic model at 
hand and can be extracted from the ground-state energy 
or 'simple' response functions Thereafter, corre- 
lation functions at long length scales or thermodynamic 
quantities at low energies can be obtained by plugging 
in K and v into analytic expressions derived within the 
exactly solvable TL model [ij, i^] . 

The recent progress in experimentally controlling iso- 
lated many-body states, in particular in cold atomic gases 
0, led to numerous theoretical studies on the dynamics 
of closed quantum systems resulting out of an abrupt 
change of the amplitude U of the two-particle interac- 
tion [8[. One assumes that the system is prepared in 
a canonical thermal state or, at temperature T = on 
which we focus, the ground state of an initial Hamilto- 
nian Hi with a given Ui; often J/j = is considered. At 
time t = the interaction is quenched to U{ and the 
time evolution is performed with the final Hamiltonian 
Hf. Fundamental questions discussed are Q: (a) Do 



some observables become stationary at large times? (b) 
How can they be classified (locality)? (c) Is it possible 
to compute their steady-state expectation values using 
an appropriate density matrix pst? It was conjectured 
that in models with many integrals of motion, e.g. those 
which are solvable by Bethe ansatz ('integrable') [9|-11|, 
Pst is not of thermal but rather of generalized ('Gibbs') 
canonical form [l2|- This has been confirmed for mod- 
els which can be mapped to effective noninteracting ones 
12" 111, in particular the TL model [l^ and a variant of 
the latter [20[ , but a general prove is lacking. For generic 
models one generally expects pst to be thermal. 

Here we address questions about the relaxation dynam- 
ics towards the steady state. Considering several models 
which in equilibrium fall into the LL class we ask: (A) 
Can the time evolution be characterized as universal in 
the LL sense? As model-dependent high energy processes 
matter at short to intermediate times one can only ex- 
pect to find LL universality on large time scales ~ but 
even then it is not obvious whether the notion of RG 
irrelevance can be transferred from equilibrium to the 
nonequilibrium dynamics pll. [22j. (B) If LL universality 
is found, does it hold independently of the number of in- 
tegrals of motion and thus independently of the expected 
nature of the steady state (generalized canonical versus 
thermal)? To answer those questions we proceed in steps. 
We compute the time evolution of the Z -factor jump of 
the momentum distribution function ri(fc, t) at the Fermi 
momentum fcp and the kinetic energy per length ekin(0 
for the spinless TL model with arbitrary momentum de- 
pendent interactions 32/4 (^) using bosonization P, 
ekin(0 is defined as the expectation value of the initial 
Hamiltonian H\ and thus describes how excitations die 
out. The quench is performed out of the noninteracting 
ground state. We analytically show that for any conti- 
nous 52/4 (^) the long-time dynamics are given by 



z{t) - 1--^^'^^) , 

|deki„(t)/dt| -e(A;v)i-3 
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The decay of Z{t) is governed by a universal exponent 
7st(i^) that depends on the equiUbrium LL parameter K 
only; ekin(i) features an asymptotic power law with an 
interaction-independent exponent but universal prefactor 
e{K, v) determined by K as well as by the renormalized 
velocity v. In equilibrium, this is the characteristic K- 
and w-dependence of correlation functions or of thermo- 
dynamic quantities, which a posteriori motivates to con- 
sider Z{t) and ekin(i) as representative examples. The 
notion of LL universality can now be defined in anal- 
ogy to equilibrium: The quench dynamics is universal if 
Eqs. ([Ip describe the long-time relaxation for any model 
falling into the equilibrium LL universality class if the 
corresponding values for K and v are plugged in. To in- 
vestigate this we compute Z{t) for a Id lattice of spinless 
fermions with nearest-neighbor hopping and interaction 
A [2^ as well as an extension of the latter including 
a next-to-nearest-neighbor interaction A2. We use the 
numerical time-dependent density-matrix RG (DMRG) 
27 1 . The model with A2 = has many conserved 
quantities, is Bethe ansatz integrable, and thus K as 
well as V are known analytically [HI, The A-A2~ 

model, however, is believed to be not exactly solvable. 
For A2 > wc extract K and v from equilibrium quanti- 
ties (e.g. the small momentum density response function) 
using DMRG [i,!!!,!!^]. Our data for the Z-factor agrees 
with Eq. ([Ta|) for any interaction strength, filling factor, 
and irrespective of the integrability of the model. The 
results for ekin(i) are consistent with Eq. (|lbp . but on 
the time scales accessible by DMRG the asymptotic be- 
havior is still masked by oscillatory terms of higher order 
in t~^. To unambiguously determine the prefactor of the 
t~'^-dccay of the energy we resort to a numerical trick. In- 
stead of performing the time evolution with exp (— i_ffft) 
we apply the imaginary time analogue exp(— TJfr). In 
this case the total energy per length e(r) - which is no 
longer conserved - is the natural observable. For the TL 
model we show that the asymptotics is completely anal- 
ogous to Eq. (fTb|) with t — r and e{K, v) replaced by 
a different function eit{K,v). For the lattice model, the 
r~^-decay manifests over several orders of magnitude, 
and the prefactor agrees with the TL prediction. 

This altogether provides strong evidence that questions 
(A) and (B) can be answered by 'yes'. We conjecture 
that the universality of the quench dynamics also holds 
for other models falling into the equilibrium LL class. 

The TL model — After bosonizing [l|, Q the density of 
left and right moving fermions with a linear dispersion 
the Hamiltonian of the TL model is quadratic in opera- 
tors bn'^ which obey bosonic commutation relations: 



n>0 



kn VF + 



27r 



+kj-^{blbl^+b_^b^) 



(2) 
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FIG. 1: (Color online) Time evolution of the ^-factor out of 
the noninteracting ground state of a Id metallic Fermi system 
after switching on two-particle terms at time t — 0. Dashed 
lines show the universal asymptotic power law t with 
an exponent determined by the equilibrium LL parameter K. 
(a) TL model. The plots displays box-like two-particle inter- 
actions g{k) of strength g = g{0); the asymptotics are univer- 
sal for any g{k). Time is given in units of (npfcc)"^. (b, c. 
Inset) Spinless lattice fermions of Eq. ^ at filling i/ featuring 
nearest (A) and next-nearest (A2) neighbor interactions. 



where kn = 2Tm/L, n ^ Tj, L denotes the chain length, 
and v-p is the Fermi velocity. The two coupling functions 
(potentials) 32/4 determine the strength of the scatter- 
ing of fermions on different branches (52) and the same 
branch (54). Usually the /c-dependence of 32/4 is ne- 
glected and integrals are regularized in the ultraviolet 
through an ad hoc procedure [H, Q • As the momentum 
dependence is RG irrelevant this is justified in equilib- 
rium if all energy scales are sent to zero Q. For the 
quench dynamics - even at asymptotic times - it is, how- 
ever, not clear if the same reasoning holds and we thus 
keep the full /c-dependence and consider coupling func- 
tions. In fact, it was recently shown that the momen- 
tum dependence indeed affects the long-time dynamics 
of certain observables 



22|. For the system to be a LL 



in equilibrium we require that < 52/4(0) < 00 (re- 
pulsive interactions) and that 52/4 (fc) decay on a scale 
kc. The Hamiltonian of Eq. ^ can be diagonalized to 
H = J2n^Q ^{kn) ctliO-n + -E-gs by introducing new modes 



= c(fc„)6„ + s(fc„)6L„ with 



i + g4(fc) 

W{k) 



\k)-l 



(3) 



Lo{k) = vp\k\ W{k) = vp\k\^{l + g^{k)Y~gl{k) , 
where 52/4 = 32/4/(27rtiF), and i?gs denoting the ground 
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FIG. 2: (Color online) Time evolution of the kinetic energy per length dckin/rfi- (a) TL model for different (Gaussian and 
quartic) two-particle potentials g{k). For any continous g{k), de^^i^/dt asymptotically falls off as e{K,v)/t'^ with a universal 
prefactor. ekin and t are given in units of vpk'^ and {vFkc)~^ , respectively, (b) Spinless lattice fermions. Solid lines show DMRG 
data; dashed lines display t^^ power laws where the ratio of prefactors is chosen according to the TL prediction. Inset: DMRG 
data before taking the t-derivative. 



state energy. The LL parameter and the renormalized 
velocity read 



y 1 + .94(0) + .92(0) 

As our initial state we take the noninteracting ground 
state I -Egg) which is given by the vacuum |vac(6)) with 
respect to the fo„. Expectation values of the time- 
evolved state \^{t)) = exp{-iHt) \E^^) can be computed 
straightforwardly using the simple time dependence of 
the eigcnmode operators and their linear dependence 
on the bn"^ 
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After bosonizing the fermionic field operator [J Q the 
Z-factor Z{t) = liiRk/^kF nik^t) — limfcx.fc^ n(k, t) is easily 
obtained (taking L oo) O [sij : 



Z{t) = exp 



/°°„ 4s2(fc)c2(fc) , r , . ,^ 

- / dk I (1 -cos[2a;(fc)^]) 

k 



Independent of the form of 52/4 (^) (even for potentials 
with a discontinous jump to zero at fee) the large-time be- 
havior is given by Eq. ((Ta]) with 7st = {K^ + K^^ — 2) /4; 
it manifests on the (nonuniversal) scale {v-pkc)~^- Fig- 
ure [TJa) shows Z(t) obtained by numerically performing 
the integral for a simple box shaped potential 52 (fc) = 
diik) ~ gQ{kc — \k\)/2 of varying amplitude g. The 
asymptotic power-law is modulated by oscillations which 
decay faster than . 

The kinetic energy per length ekin(i) reads (L 00) 



ekin(i) 



Vp 



dkkAs^{k)c^{k) {1 - cos [2uj{k)t]}. (5) 



The steady-state value is obtained by dropping the os- 
cillatory term which averages out for i — > 00. For con- 
tinuous coupling functions 52/4 (fc) of range kc asymp- 
totic analysis yields Eq. (jlb[) as the leading term in the 



long-time limit: j33j the coefficient is given by e(K,v) = 
7st(i^)wF/(47ru^). Figure m^a) shows the derivative of 
ekin for 52(A) = 54(fc) = g{k), a Gaussian potential 
g{k) = gexp(— [fc/fcc]^/2)/2 as well as a quartic poten- 
tial g{k) = .9/(1 + [k/kc]'^)/2 and varying interaction 
strengths. As either g{0) or the lowest nonvanishing Tay- 
lor expansion order of g{k)—g{0) increases, the amplitude 
of an oscillatory term which decays faster than the lead- 
ing one becomes stronger. The (nonuniversal) scale on 
which the asymptotic t~^-behavior dominates thus heav- 
ily depends on the strength and type of potential at hand 
[compare the inset and the main part of Figure [^a)]. 

Microscopic lattice model — As a next step we provide 
strong evidence that Eqs. ([T|) describe the long-time re- 
laxation dynamics of any model which in equilibrium falls 
into the LL universality class. To this end, we consider 
spinless lattice fermions. 



H- 



1 1 



(6) 



with Uj = c^jCj — 1/2. We study the quench dynamics 
using an infinite-system DMRG algorithm 25|, |32| . We 
determine l^'gs) by applying an imaginary time evolution 
exp(— r7?|A=A2=o) to a random initial matrix product 
state with a fixed matrix dimension x until the energy 
has converged to typically 8 — 10 relative digits. Oper- 
ators exp(~ H) are factorized by a second or fourth or- 
der Trotter decomposition. Thereafter, we compute the 
real time evolution |^'(i)) = exp{—itH)\E^g) in presence 
of the two-particle terms A and A2. x is dynamically 
increased in order to maintain a fixed discarded weight. 
We carefully ensure that the latter is chosen small enough 
(and that the initial x is large enough) to obtain numer- 
ically exact results. 

The time evolution of the momentum distribution 
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(Color 



online) Imaginary time evolution 
exp(— rff)|_Ega) ^-^4° l-Egs) towards the interacting ground 
state l-Egs). The main part shows DMRG data for the 
T-derivative of the total energy in the lattice model. The 
T~'^-decay predicted by bosonization manifests over several 
orders of magnitude. The prefactor agrees with the TL 
formula ^\x(K,v) for all parameters (this is illustrated in 
the insets). The imaginary-time energy dynamics are thus 
universal in the LL sense. 



function n(fc,i) = Ej e'''-'(*W|c]sl*W) 
responding Z-factor can be computed straightforvifardly; 
the j-sum is carried out up to ^ 10000 sites. Results for 
Z(t) are shown in Figure [ljb,c). Its long-time asymp- 
totics indeed shows a power-law decay and the 

exponent agrees to the one predicted by the TL formula 
(dashed lines) if the latter is evaluated for K correspond- 
ing to the microscopic parameters under consideration. 
We take K from the Bethe ansatz ( A2 = 0)0, [11] or the 
equilibrium density response (A2 > 0) [l,!!!]- The agree- 
ment with the TL model result holds for any interaction 
strength (soj . for any filling factor and irrespective of 
the integrability of the model. This strongly indicates 
that the asymptotic dynamics of the Z-factor is indeed 
universal in the LL sense. 

The time derivative of the kinetic energy (the expec- 
tation value of -ff |a=A2=o) per length (site) is shown in 
Figure [2jb). Its magnitude decays as ■, and the ra- 
tio between prefactors (where the factor drops out) 
at different interaction strengths is consistent with the 
TL formula (see the dashed lines). However, oscilla- 
tions have not died out completely, and a pure power 
law cannot be identified unambigously. To further sup- 
port that this is merely because the time scales reach- 
able in our DMRG calculation are too small - remind 
that for the TL model the scale where deYmjAt is gov- 
erned by Eq. (jlbp strongly depends on the strength and 
type of the potential in contrast to Zify where it is al- 
ways (upfcc)"^ - and that the energy relaxation is in- 
deed universal, we consider an imaginary time evolution, 

\nr)) = exp(-7Jr)|i?o,)/(i?o,|exp(-2ifr)|£;0,)^/^ 



of the interacting Hamiltonian. The total energy is the 
natural observable to compute in this academic scenario. 
Its asymptotic behavior within the TL model is com- 
pletely analogous to Eq. (jlbp with i — > t and e{K, v) 
eit{K,v) = U2{[K + K-^ - 2]/[K + + 2])/(8^w), 
where Li2 denotes the dilogarithm [s^l- For the lattice 
model, one can easily access large imaginary times using 
DMRG, and the T~^-decay manifests over several orders 
of magnitude. This is illustrated in Fig. [3l The prefac- 
tor (shown in the insets) agrees with en {K, v) (the latter 
depends on K and v only; thus, one does not need to 
consider ratios) for all interactions and fillings. The dy- 
namics of the total energy at large r is universal. 

Conclusion — We have obtained exact expressions for 
the time evolution of the Z-factor and of the kinetic en- 
ergy after an interaction quench within the Tomonaga- 
Luttinger model. For any continous two-particle poten- 
tial, their long-time asymptotes Z ^ t^'''^', dekin/c^^ ~ 
e{K, v)/t^ are universal functions of the LL parameter K 
and the renormalized velocity v. We studied a similar 
scenario for spinless lattice fermions using DMRG; for 
large times, Z{t) and ekin(0 are described by the above 
expressions. This provides strong evidence that the re- 
laxation dynamics after an interaction quench within any 
model that falls into the equilibrium LL universality class 
has aspects which are universal in the LL sense. 
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